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Abstract 

Recently, there has been growing interest in employing condensed matter systems such as quan- 
tum spin or harmonic chains as quantum channels for short distance communication. Many prop- 
erties of such chains are determined by the spectral gap between their ground and excited states. 
In particular this gap vanishes at critical points of quantum phase transitions. In this article we 
study the relation between the transfer speed and quality of such a system and the size of its 
spectral gap. We find that the transfer is almost perfect but slow for large spectral gaps and fast 
but rather inefficient for small gaps. 
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I. INTRODUCTION 



Quantum Information and Condensed Matter Physics are currently two areas of very 
active research where several links between the two branches have been discovered in recent 
years. In particular, the relation of entanglement properties to critical points of quantum 
phase transitions has been studied in some detail. In analogy to classical phase transitions, 
the latter are typically characterized in terms of the scaling behavior of equilibrium proper- 
ties, such as the correlation length. Analogous scaling phenomena have now been found for 



m. 



the entanglement properties of a spin chain in the vicinity of a critical point. 

In Quantum Communication, condensed matter systems have recently received some 
attention as interesting candidates for quantum channels . This approach may be of 
advantage in situations where photonic quantum communication is not possible for example 
because the two units are only separated by a few optical wavelengths or are inside a material 
that does not permit for the well-controlled propagation of light. 

The quantum mechanical properties of correlated many body system are to some extent 
determined by the energetic gap between their ground states and the lowest excited states. 
In particular the system is short range correlated if it features a finite gap, while some 
quantities can become long range correlated when the gap closes 0,Q|. The latter happens 
at the critical points of quantum phase transitions jlj]. 

Motivated by these findings the dynamical entanglement properties of quantum many 
body systems undergoing a quantum phase transition are receiving increasing attention 
10, U, 12, 3- 

To pursue the analysis of quantum many body systems as possible quantum channels 
a step further, we study here the connection between the spectral gap and the systems 
ability to transfer excitations and entanglement. However, our aim is twofold: Besides 
being interested in the systems transfer capacity, we also want to explore, to what degree 
the transfer properties can be employed to detect and characterize critical points of quantum 
phase transitions. 

In particular, we consider spin chains and harmonic chains, where in both cases, two 
ancillas couple weakly to the chain at distant sites. We study the transfer of quantum 
information and excitations from one ancilla to the other numerically, by employing newly 
developed matrix product state techniques Q|, and analytically with a master equation 
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approach. The numerics is applied to the spin chains while the master equation is used 
to describe the case of an harmonic chain, where its validity is confirmed numerically for 
Gaussian initial states Q|. 

Both systems show a very similar behavior: The transfer properties sensitively depend 
on the energy gap between the ground state and the lowest excited states, but do not sig- 
nificantly dependent on the detailed structure of the Hamiltonian. The quality and speed of 
state transfer through such systems may thus be used to detect critical points experimentally. 

This paper is organized as follows: Section II begins with a discussion of spin chains, 
mainly employing numerical methods. In Section III a simplified picture unraveling the 
origin of the observed behavior is presented. This picture is then further corroborated 
in Section IV where we study harmonic chains both numerically and analytically with a 
master equation approach. Section V discusses connection to quantum channel capacities 
and Section VI summarizes our results. 



II. SPIN CHAINS 



A variety of quantum spin chains exhibit quantum phase transitions. Therefore we be- 
gin with the study of a linear chain of spins with nearest neighbor interactions and open 
boundary conditions described by the Hamiltonian 

N N-l 
i=l i=l 

Here N is the number of spins, B a magnetic field and J x , J y and J z the interaction between 
neighboring spins. Two ancilla systems (named S for "sender" and R for "receiver" ) couple 
to the chain at spins ms and mjj, which are near the center of the chain in order to avoid 
boundary effects. The total Hamiltonian of chain and ancillas then reads 

H = i/ cham + B a K + o& + J a (o*goZ B + o%pZ H ) , (2) 

where B a > is the Zeeman splitting of the ancillas, which might differ from B, and J a > 
the coupling of the ancillas to the chain, which is taken to be weak, i.e. J a (B, J x , J y , J z ). 
As an initial state of the system we consider 

|*(0)) = |Ts4fl,0cham>, (3) 
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FIG. 1: P(lslR)(t) (light gray), P(UU)(t) (gray) and P(isU)(t) (black) for B = 1, J x = 0.3, 
J, -J.- 0, B. - 0,4 and J. - 0.05 as given by th. —ion for th. open boundary .nodel 
with N = 100 spins. 5 couples to spin 45 and R to spin 55. Figure taken from [131 ] . 



i.e. the sender is spin up and the receiver is spin down, while the chain is in its ground state, 
|0 c hain)- This will allow us to explore the propagation of this excitation through the chain. 

The dynamics is simulated with a matrix product state technique using matrices of dimen- 
sion io x io QQ. The accuracy of the simulations was confirmed by varying the matrix 
dimension and the size of the time steps. Furthermore, we confirmed that the dynamics 
conserves the total energy of the system. 

Figure ^ shows the probability -P(TsIb) that "sender" S is in its excited state | \g) and 
the "receiver" R in its ground state | together with P(Is1r) and P(IsIr) for a model 
with N = 100, m s = 45, m R = 55, B = 1, J x = 0.3, J y = J z = 0, B a = 0.64 and J a = 0.05. 
-PQsT-r) is always less than 10 -4 . The plots show that the excitation that was initially 
located in S oscillates back and forth between S and R. 

Figure shows P(tsU), P(isU) and P(IsIr) for a model with N = 600, m s = 295, 
m R = 305, B = 1, J x = 0.3, J y = J z = 0, B a = 0.8 and J a = 0.05. Again, P(tsTfl) 
is always less than 10 -4 . For these parameters, the excitation is not fully transferred to 
R, contrary to figure ^ Both, S and R relax to their ground states with the excitation 
only being partially and temporarily transferred to R, even for close-lying spins. Note that 
the parameters chosen in figures Q and El are the same except for B a which in figure El is 
significantly larger than in figure ^ 

The two observed scenarios are rather generic. To demonstrate this, we have done the 
same simulations for a XXZ-model. The results, shown in figures El and 01 clearly agree with 
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FIG. 2: P(lslR)(t) (dotted line), P(Ts4r)W (dashed line) and P(lsU)(t) (solid line) for B = 1, 
J x = 0.3, J y = J z = 0, B a = 0.8 and J a = 0.05 as given by the simulation for the open boundary 



model with iV = 600 spins. S couples to spin 295 and R to spin 305. Figure taken from 
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FIG. 3: P(UU)(t) (light gray), P(Mil)(<) (gray) and P(| S j H )(t) (black) for 5 = 1, J x = 0.5, 
, s . 0,, . 0,, B„ . 0.04 and . 0.5 a, given by the simU la ti o„ f„ r th. open bounds 
model with iV = 100 spins. 5 couples to spin 45 and R to spin 55. Figure taken from t 13]. 

our findings for the previous coupling parameters. Again B a in figure 0] is significantly larger 
than in figure EJ while all other parameters are equal. 

III. SIMPLIFIED PHYSICAL PICTURE 



The exact dynamics of a quantum spin chain is extraordinarily complex. Nevertheless, 
underlying the dramatic difference between the almost perfect transfer scenarios in figures 
Q and El and the damped scenario in figures El and |U is a simple physical mechanism. The 
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FIG. 4: P(isi R )(t) (dotted line), P(Tslfl)(*) (dashed line) and P(lsU){t) (solid line) for B = 1, 
Jr = 0.3, J y = 0.2, J z = 0.1, -B a = 0.2 and J a = 0.05 as given by the simulation for the open 
boundary model with N = 600 spins. S couples to spin 295 and R to spin 305. Figure taken from 

dynamics we have simulated is given by the Schrodinger equation containing the Hamiltonian 
(J2J). As a consequence, all moments of the Hamiltonian are conserved. The initial state 
|\l/(0)), where one of the ancillas is excited, is not an eigenstate of H as given by The 
mean value is (\I / (0)|iJ|\I/(0)) and the variance of the energy is given by 



Given that the evolution is assumed to be Hamiltonian, both are in fact time invariant. Thus 
in the entire dynamics only those states with an energy expectation value E in the range 
(*(0)|#|*(0)) - V2J a <E< (*(0)|fT|*(0)) + V2J a play a significant role as indicated 
by Figure 03 S and R are depicted as two level systems, while for the chain there is a 
unique ground state and a quasi continuous band of excited states sketched as the gray 
area. The dots indicate the initial occupations. The relevant energy range given by the 
variance lies between the two horizontal dashed lines. If the spectral gap is larger than the 
Zeeman splitting of the ancillas (left plot) it is not possible to create a real excitation of the 
chain, only processes involving virtual excitations of the chain are then permitted and the 
excitation will be almost completely transferred from S to R. If however the spectral gap 
is smaller than the Zeeman splitting of the ancillas, real excitations may be created in the 
chain. These excitations dephase rapidly due to the large number of accessible energy levels 
which prevents their reabsorption - the chain acts as a bath and information is lost. 



AE = (0)1 H 2 |tt(0)> - (^(0)|#|^(0)> 2 = V2J a . 



(4) 
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FIG. 5: Sketch of the energy levels of the system. The dots indicate the occupations of the initial 
state. For this initial state only the energy levels between the two horizontal dashed lines are 
accessible, resulting in almost perfect transfer for the left scenario and damping for the right one. 

We believe that this heuristic picture captures the essential physics. To further corrobo- 
rate our simple physical model and to underline the generality of our findings, we now turn 
to a different model for the chain which also features an adjustable energy gap above its 
unique ground state. 

IV. HARMONIC CHAIN 

Harmonic chains whose time evolution is governed by a Hamiltonian that is quadratic 
in position and momentum can be solved exactly in a compact form and therefore lend 
themselves to analytical investigations. Furthermore, Gaussian initial states, ie states whose 
Wigner function is a Gaussian, preserve their Gaussian character under the time evolution. 
As a Gaussian is fully described by its first and second moments this allows for a very 
compact, polynomial in the number of oscillators, description of the state of the systems 
and its evolution in time (see e.g. for details), which permits the numerical investigation 
of very large systems. This is the basic motivation for the study of harmonic systems in 
Gaussian states in this section. 

We consider a harmonic chain with periodic boundary conditions described by 



with the Pj being the momenta and the qj the positions (qN+i = qi). In this case, the two 
ancillas are harmonic oscillators that couple to oscillators ms and m# of the chain. The 





complete Hamiltonian reads 



Htot — -f^chain + -^ancillas + Hj (6) 
^ancillas = ~ (p| + + P R + ^Qr) (7) 

Hi = JaiqsQms + <?£!<?mj • (8) 

We are interested in the time evolution of the ancillas only. Thus, we will now derive an 
approximate master equation for the reduced density matrix a(t) of the two ancilla systems. 
For weak coupling J a <C (Q, fl ) we find 



da 
~dt 



[ dsTr chain {[H I (t),[H I (s),\0)(0\®a(s)}}} , (9) 
J o 



where a(t) is the reduced density matrix of the ancillas, Tr chain is the trace over the degrees 
of freedom of the chain and |0) denotes the ground state of the chain (|5|). Hj(t) is the 
interaction between ancillas and chain in the interaction picture. Here H = Hq + Hi with 

Hq = -ffchain + -^ancillas and 

H^t) = e iHat H ie - iHot 

and 

a{t) = e iHot pe- iHot . 

Eq. © is an expansion in the coupling strength J a up to second order, which is a good 
approximation if the integral approaches a constant value for t > t*, where J a t* <C 1. Since 
a only changes significantly on time scales t ~ J^ 1 ^> t*, the approximation a(s) m a(t) 
can be used. Performing the trace on the rhs of Q yields 



da 
~dt 



-Ja E (-i(Yx + (Y -Y 1 )6 jl )[a j ala] 

+ (X x + (X - X 1 )S jl ) ({a ja ta} - 2{ aj aa\))) , (10) 



j,l=S,R 



where [•, •] and {•, •} denote commutators and anti-commutators and as and o,r are the 
annihilation operators of S and R, respectively 



1 

2uj 
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for j = S,R. In deriving (fTUj) . we applied a rotating wave approximation. The coefficients 
read 



1 
1 



2u 



^1 2 ( ^ m ^ ,rls?ni? ^ m S™R' 



C£(*)= / <M0|ft(0a(*)|0> (11) 



with given by 



where fc, / = m s ,m R . Due to the symmetries of the model, we have Cj^ sms = C^ flmfl 
and Cm smfl = C^ lRms . Eq. (fT0|) is a good approximation whenever Cj^(i) « C w = 
const, for t <C J^ 1 . Since the C kl (t) do not depend on J a themselves, there is always 
a sufficiently small J a such that this holds, provided lim^oo C kl (t) exists. The validity of 
the master equation can be confirmed by a numerical simulation for a chain with 1400 oscil- 
lators. We found good agreement between our analytical and numerical solutions, with the 
relative errors being less than 5%. 

The _ chain can be diagonal^ via a Fonrier transfer Q In the linnt of an 
infinitely long chain, N — > oo, the correlation functions read 

(0|^(*) ft (a)|0> = ^ f dk— oos((; -l)k)e-^- s \ (12) 

^ JO 

where <jj\ = 4f2 2 sin 2 (/c/2) + Q 2 , with — 7r < k < tt. and all lim^oo C^(t) exist except for the 
case where u; = Qq = 0. As in master equations for system bath models, we now insert the 
asymptotic expressions, C% = lim^oo C kl (t), into eq. (fTUj) . This replacement assumes that 
all internal dynamics of the chain happens on much shorter time scales than the dynamics 
caused by the interaction of the ancillas with the chain. Furthermore, it does not treat the 
initial evolution for short times with full accuracy since lim^o C^(t) = C kl ). The master 
equation thus obtained is therefore valid in a regime where the couplings J a are weak enough 
such that the time it takes for an excitation to travel from S to R is completely determined 
by J a , i.e. by the time it takes to be transferred into and from the chain. 

From eq. (fTU|) we find the following solution for the expectation values of the occupation 
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FIG. 6: Frequencies 2J^j/i of the excitation's oscillation between S and R for $7 = 1, J a = 0.05 



and uj = 0.5. The transfer speed increases as — w — > 0. Figure taken from 
numbers of 5* and R, ns = Tr (a^ascr) and n R = Tr(a^ORcr): 



M0 1 , (M0) + M0) cosh{2j!xit) ± MQ)-^(Q) cos(2 ^ A e - 2 ^« (13) 



Here, a; = lim^oo X , X\ = lim^oo X\ and y\ = lim^oo Y 1 . Note that x > and x > \x\\. 
Inserting the correlation functions into eq. (fTTj) one sees that Xq and x\ are only non-zero if 
it) > Uk for at least one mode k, that is if our initial state is in resonance with states where 
both ancillas are in their ground states and the chain is in one of its lowest-lying excited 
states (c.f. figure El). The dispersion relation shows that this only happens for uj > Q . As 
in figures ^ and 121 or El and ^ we thus observe two different scenarios: 

If uj < Qq, and therefore Xq = X\ = 0, the excitation that is initially in S oscillates back 
and forth between S and R at a frequency 2J^yi, i.e. 

n s {0)+n R (0) ,n s {0)-n R {0) fn l2 




± " V \ UK ' cos(2J Q V)- (14) 



Note in particular that the excitation is entirely transferred to R at times t n = 
n{it j JaVi)'i n — 1,2, ... . Figure El shows the frequencies 2J^yi of the excitation's oscil- 
lations between S and R for cases where uj < Qq, for Q = 1, J a = 0.05 and uj = 0.35 
as a function of Qq. As Qq — uj decreases, the transfer becomes faster and the oscillation 
frequency increases. 

Let us examine the dependence of the oscillation frequency 2J„yi on the distance \ms — 
m R \. With the definitions of y\ and Y 1 and eqs. (fTTf . ([T2j) and J °° dr e~ lXT = n5(0) — %P\ , 
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where V denotes the principal value, we find 

Using the substitution z = e lfc , this integral can be converted to an integral over the unit 



circle, which in turn can be evaluated via the residue theorem 1151. The result is 



j2 z \m s -m R \ _|_ z ~\ms— w,r\ 



2J ^ = -sfe — to — (16 » 



where z = 2y/a 2 + a — (2a + 1) and a = (Qq — u 2 )/4Q 2 . The frequency 2J 2 yi and hence 
the transfer speed decrease exponentially with increasing |m$ — m^|. 

If, on the other hand, u > f2o, the chain acts similarly to a bath. Here, Xo 7^ 0, X\ 7^ 0, and 
both ancillas relax into their ground state transferring their energy into the chain. During 
this process, however, a fraction of the energy initially located in S appears momentarily in R 
before it is finally damped into the chain. The maximal excitation of the receiver throughout 
the entire evolution depends on the distance \ms — mn\. For a given initial energy in S, only 
a narrow range of the excitation spectrum is relevant for the dynamics. Hence all relevant 
modes of the chain have very similar wavelength. If the distance \ms — turI is a multiple of 
half that wavelength, a wave that has high amplitude next to S will have high amplitude 
next to R, too. In that case, the maximal excitation probability in R will be relatively large, 
although less than 0.5. If on the other hand, a wave that has high amplitude next to S 
has zero amplitude next to R, the transfer vanishes completely. Figure [7| shows the solution 
(THH) for a harmonic chain and ancillas with Q = 1, u = 0.5, J a = 0.05, |mg — = 9 and 
Qo = 0.2. One might try to derive the same type of master equation for the spin chain (pj. 
In the next section, we explain why such an approach is problematic when attempted with 
standard techniques. 



A. The problem of even and odd excitation number subspaces for spin chains 

The standard procedure to diagonalize say the transverse Ising chain is a Jordan Wigner 
transformation followed by Fourier and Bogoliubov transformations ^j. This procedure 
however only diagonalizes the subspace of odd or even number of excitations at a time but 
not both. Since the operator a? maps between the odd to the even number subspaces for any 
j, an exact expression for the time dependent correlation functions (0 1 crj (t) <rf (s) 1 0) , which 
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FIG. 7: The solution (1131) for a harmonic chain and ancillas with Q, = 1, u) 
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\ms — mn\ = 9 an d = 0.2. Figure taken from |18l |. 
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FIG. 8: P(\s[R){t) (dashed line) and P(|stfl)(*) (solid line) for B = 1, J x = 0.3, J y = J z = 0.0, 
B a = 0.8 and J a = 0.05 as given by the solution of a master equation for the transverse Ising spin 
chain. The solution is not physical since probabilities are not less or equal to 1. 

would be required for deriving a master equation, cannot be obtained via the above procedure 
[bs| . If one uses only say the even excitation number subspace as an approximation, the 
errors for the obtained eigenvalues scale as 1/N, where N is the length of the chain. It is 
therefore often argued that only using one subspace is a good approximation for long chains. 

In our problem, this reasoning cannot be applied: An approximation of the required 
correlation functions, (0\(7j(t)of(s)\0), by only considering one subspace leads to unphysical 
results for the master equation. An example is given in figure El We have therefore avoided 
a treatment of the spin chains with master equations. 
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V. CONCLUSIONS 



Quantum spin chains may be used as quantum channels. Here we have studied the 
dependence of the transfer quality and speed on the size of the spectral gap between the 
ground and the lowest excited state of the considered system. 

Two different scenarios are observed: Provided the spectral gap of the chain is larger than 
the amount of energy initially located in the sender, close to perfect transfer happens on 
rather long time scales. The transfer becomes even more perfect but slower if the coupling 
strength between ancillas and chain or the energy initially in the sender is lowered further. 

If the energy in the sender is initially larger than the spectral gap, the behavior changes 
significantly. Here, the excitations get lost into the chain implying a bad transfer quality. 
This second scenario takes place on much shorter time scales than the first one. 

These results may be used for two main applications: For quantum information purposes, 
they show, that by using a system with a large enough spectral gap one can build good 
quantum channels. These channels are furthermore considerably robust against decoherence 
because the chain remains to a high precision in its ground state which is not significantly 
affected by an environment at zero temperature. 

In condensed matter physics, on the other hand, the transfer speed and quality can be 
used to infer the size of the spectral gap and hence the location of quantum critical points. 
This method is admittedly inferior to neutron scattering in magnetic composites, but it may 
by advantageous in optical lattices where neutron scattering is not that successfull. 

Acknowledgments 

The authors would like to thank Sougato Bose and Daniel Burgarth for discussions at 
early stages of this project as well as comments by Andrew Fisher. This work is part of the 
QIP-IRC supported by EPSRC (GR/S82176/0), the Integrated Project Qubit Applications 
(QAP) supported by the 1ST directorate as Contract Number 015848' and was supported 



13 



also by the Alexander von Humboldt Foundation, Hewlett-Packard and the Royal Society. 



[1] Bose S. Phys. Rev. Lett. 91, 207901 (2003). 

Burgarth D and Bose S. New J. Phys. 7, 135 2005. 
[2] Plenio M B, Hartley J, and Eisert J. New J. Phys. 6, 36 (2004). 

Eisert J, Plenio M B, Bose S, and Hartley J. Phys. Rev. Lett. 93, 190402 (2004). 

Plenio M B and Semiao F. New J. Phys. 7, 73 (2005). 

Perales A and Plenio M B J. Opt. B 7, 601 (2005) 
[3] Hastings M B Phys. Rev. Lett. 93, 126402 (2004). 

Hastings M B Phys. Rev. B 69, 104431 (2004). 

Hastings M B Phys. Rev. B 73, 085115 (2006). 

Hastings M B Phys. Rev. Lett. 93, 140402 (2004). 
[4] Cramer M, Eisert J, Plenio M B and Dreissig J. Phys. Rev. A 73, 012309 (2006). 

Cramer M and Eisert J. New J. Phys. 8, 71 (2006). 

Wolf M M, Giedke G and Cirac J I. Phys. Rev. Lett. 96, 080502 (2006). 
[5] Sachdev S Quantum Phase Transitions. Cambridge Univ. Press, Cambridge, (1999). 
[6] Osterloh A, Amico L, Falci G and Fazio R. Nature 416, 608 (2002). 
[7] Verstraete F, Martm-Delgado M A and Cirac J I. Phys. Rev. Lett. 92, 087201 (2004). 

Pachos J K and Plenio M B. Phys. Rev. Lett. 93, 056402 (2004). 

Kay A, Lee D K K, Pachos J K, Plenio M B, Reuter M E and Rico E. Optics and Spectroscopy 
99, 339 - 356 (2005). 
[8] Larsson D and Johannesson H. Phys. Rev. Lett. 95, 196406 (2005). 

[9] Amico L, Osterloh A, Plastina F, Fazio R and Massimo Palma G Phys. Rev. A 69, 022304 

(2004) . 

Diir W, Hartmann L, Hein M, Lewenstein M and Briegel H-J Phys. Rev. Lett. 94, 097203 

(2005) . 

Sen De A, Sen U, Ahufinger V, Briegel H-J, Sanpera A and Lewenstein M quant-ph/0507172 
(2005). 

[10] Sen A, Sen U and Lewenstein M. Phys. Rev. A 72, 052319 (2005). 
[11] Yi X X, Cui H T and Wang L C. |quant-ph/0511026| (2005). 



11 



[12] Zurek W H, Dorner U and Zoller P. Phys. Rev. Lett. 95, 105701 (2005). 

Polkovnikov A. Phys. Rev. B 72, 161201 (R) (2005). 

Dziarmaga J. Phys. Rev. Lett. 95, 245701 (2005) 
[13] Hartmann M J, Reuter M E, and Plenio M B. New J. Phys. 8, 94 (2006). 
[14] Ostlund S, Rommer S Phys. Rev. Lett. 75, 3537-3540 (1995). 

Vidal G. Phys. Rev. Lett. 93, 040502 (2004). 

Daley A J, Kollath C, Schollwock U and Vidal G. J. Stat. Mech.: Theor. Exp., P04005 (2004). 
[15] We like to thank Andrew Fisher for bringing this point to our attention. 
[16] Plenio M B, Eisert J, Dreissig J, and Cramer M. Phys. Rev. Lett. 94, 060503 (2005). 
[17] Lieb E, Schultz T and Mattis D. Ann. Phys. 16, 407 (1961). 
[18] McCoy B M. Phys. Rev. A 4, 2331 (1971). 
[19] Eisert J. and Plenio M.B., Int. J. Quant. Inform. 1, 479 (2003) 



15 



